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O ■ ABSTRACT 

(N 

^ . We consider hot accretion onto a rapidly spinning neutron star (or any other 

^^ , compact object with a surface). A radiative hot settling flow has been discovered 

at low accretion rates in the early work by Medvedev & Narayan (2001) and 
analytical solution has been presented. It was shown later that this flow can 
match external medium smoothly, thus enforcing its physical feasibility. Here we 



> 

t^^ . complete the study of the global structure of such hot accretion by presenting the 



analytical solution for the boundary later, which forms between the bulk of the 



m 

O ! flow and the stellar surface. We confirm our results via full numerical solution of 

Q ! height-integrated two-temperature hydrodynamic equations. 

H . 1. Introduction 

^ ' Accretion flows around compact objects frequently radiate significant levels of hard X- 

rays, indicating the presence of hot optically-thin gas in these systems. This has motivated 
/\ I the study of hot accretion flows around compact stars. 

Various models of accretion may be divided into two large groups according to whether 
they require a shock to form near the central object (e.g., a neutron star) or not. The first 
group includes spherical accretion flows onto a neutron star (NS) and/or a black hole (BH) 
in either a kinetic (Zeldovich & Shakura 1969; Alme & Wilson 1973; Turolla et al. 1994; 
Zampieri et al. 1995; Zane et al. 1998) or a fluid (hydrodynamic) regime (Shapiro & Salpeter 
1975; Chakrabarti & Sahu 1997; Kylafis & Lamb 1982). Whereas the existence of a shock 
is not a problem in spherical flows, the shock discontinuity in rotating flows likely leads to 
serious physical inconsistencies associated with causality of viscous processes (Pringle 1977; 
Popham & Narayan 1992) [except for the case of "gap accretion" in which the surface of the 
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accreting star lies below the marginally stable orbit and, hence, is not in causal contact with 
the bulk of the flow (Kluzniak & Wilson 1991; Deufel, Dullemond & Spruit 2001)]. 

The other group includes flows which contain no shocks, which (or, at least, the inner 
regions of which) thus, must be subsonic and have a gas heated up to the nearly virial 
temperature. Such hot boundary layers form in the cool thin disk accretion, provided a causal 
viscosity prescription is used (Popham & Narayan 1992; Narayan & Popham 1993; Popham & 
Sunyaev 2001). A large class of hot accretion flows: the SLE solution (Shapiro, Lightman, & 
Eardley 1976), the advection-dominated accretion flow (ADAF) (Ichimaru 1977; Rees et al. 
1982; Narayan & Yi 1994, 1995a,b; Abramowicz et al. 1995), the advection-dominated inflow- 
outflow solution (ADIOS) (Blandford & Begelman 1999), and the convection-dominated 
accretion flow (CDAF) (Narayan et al. 2000; Quataert & Gruzinov 2000). The relevance 
of these solutions for accretion onto a NS is unclear. Moreover, since the marginally stable 
orbit appears not to play an important role in hot quasi-spherical flows (Narayan, Kato, & 
Honma 1997; Chen, Abramowicz & Lasota 1997), it is, therefore, not clear that one would 
necessarily have a shock with a hot flow even if the accreting star were very compact. Finally, 
a subsonic hot accretion flow may form around a magnetized neutron star in the propeller 
state (Davies & Pringle 1981; Ikhsanov 2001, 2003) [in which the gas heating by viscosity is 
much faster than radiative cooling] and the "hot settling flow" (Medvedev & Narayan 2001; 
Narayan & Medvedev 2003) [in which heating and cooling balance each other] deflnitely fall 
in this group. 

The latter, hot settling flow, could equally well be described as a "hot atmosphere" 
since the solution is, to flrst approximation, static, and accretion represents only a small 
perturbation on the static solution. To our knowledge, the hot settling flow is the only 
steady state solution for accretion onto a NS presently available that does not involve a 
discontinuity near the surface of the star. The hot settling flow exists at rather low accretion 
rates, M, smaller than a few percent of Eddington. The flow is subsonic everywhere; it is 
powered by the rotational energy of the central accretor which is braked by viscous torques. A 
very interesting property of the flow is that, except for the inflow velocity, all gas properties, 
such as density, temperature, angular velocity, luminosity, and angular momentum flux, are 
independent of the mass accretion rate; the flow properties are sensitive to the star spin, 
s (Medvedev & Narayan 2001). Because too few parameters control the flow structure 
(effectively, just M and s) it was not clear whether the flow can match to a general external 
medium (for a flxed spin of a star). Recently, Narayan & Medvedev (2003) demonstrated 
that the flow matches the external medium via a special "bridging solution" , which forms in 
a relatively narrow region near the outer boundary. 

To complete the study of the structure of hot settling accretion flow onto a neutron star. 
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it is, thus, required to study how this flow matches the inner boundary conditions set by the 
parameters of the star. The study of the inner boundary layer is presented in this paper. 
We present an analytical self-similar solution and conflrm it numerically. Our results are in 
agreement with other studies (Titarchuk, Lapidus & Muslimov 1998; Titarchuk & Osherovich 
1999), and may be important for the interpretation of kHz quasi-periodic oscillations. We 
critically examine the limitations of our model as well. 



2. The model 

We consider viscous hydrodynamic accretion onto a compact spinning object with a 
surface. The central object has a radius i?* = r^^Rg (where Rg = 2GM^/c'^ = 2.8 x lO^m cm 
is the Schwarzchild radius), a mass M* = ttiMq, and an angular velocity fi* = sQk*, where 
D.k{R) = {GM^I R^y/'^ is the Keplerian angular velocity at radius R and VLk* = VtK{R*)- 
The mass accretion rate is M = rhMEdd, where MEdd = 1-39 x lO^^m g s^^ (corresponding to 
a radiative efficiency of 10%). We use the height-integrated form of the viscous hydrodynamic 
equations (Ichimaru 1977; Abramowicz, et al. 1988; Paczyhski 1991; Narayan & Yi 1994). 

In standard accretion problems the radial coordinate, i?, varies through a large range. 
Hence an analytical self-similar solution, in which the gas parameters (density, temperature, 
and such) scale as power-laws of i?, is usually possible in the region far from the boundaries. 

Unlike the bulk of the flow, the solution to the boundary layer (BL) cannot be obtained 
in a self-similar form in terms of the radial coordinate R. Indeed, the structure of the BL is 
intrinsically non-self-similar in R because all the gas parameters (e.g., the temperature, gas 
density, etc.) change very dramatically over a relatively short radial region: -R* < -R < 2/?^,. 
For instance, the density nearly diverges as one gets close to the star surface whereas the 
temperature decreases to the values well below the virial temperature. Such a behavior, 
however, suggests to look for a self-similar solution in terms of the distance from the stellar 
surface, i.e., in terms of 

D = R-R,. (1) 

In our calculation we neglect the effects of radiation transfer and Comptonization. They 
may be important in hot regions, but will unlikely strongly affect the flow closer to the star, 
where the temperature of the gas falls below fewxlO^ K or so (see discussion in Section 4). 

Unlike the hot settling flow case, here we cannot neglect the radial (infall) velocity. We 
use the height-integrated two-temperature hydrodynamic equations, written in the approxi- 
mation that R = R^ + D with D <^ R^, : 

-M = AttR^pv, (2) 
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,%^iff-nl.)R,^l±(^^:), (3) 

4vraK)^^ = J-Mi?^, (4) 



where p = Pp + Pe ^ Pp is the mass density of the accreting gas (pp and pe — (^e/'"^p)p 
are the mass densities of the proton and electron fluid), v is the radial infall velocity [note 
that in equation (2) we took into account that the radial velocity is negative (inward)], Q 
is the angular velocity, cj. = cj.p + cj.^ rrie/mp is the square of the isothermal non-relativistic 
sound speed [this follows from the fact that both electrons and protons contribute to the 
total gas pressure: pc^ = p = Pp + Pe = PpCsp + Pecl^ ~ p(c^p + cl^me/mp)] and Cgp and 
Cse are the respective isothermal sound speeds of the two species, a is the Shakura-Sunyaev 
viscosity parameter, J=const. is the rate of accretion of angular momentum (an eigenvalue 
of a problem), 7p and 7e are the adiabatic indexes of protons and electrons, which we assume 
to be equal, for simplicity. We assumed that the flow is hot and quasi-spherical, i.e., the 
local vertical scale height H = Cs/^k is comparable to the local radius R c::^ Rt^. We may 
then set, for simplicity, that H = R (Medvedev & Narayan 2001). Finally, g^, g~, and gcoui 
are the viscous heating rate, radiative cooling rate, and energy transfer rate from protons to 
electrons via Coulomb collisions, per unit mass. We have assumed that all viscous heat goes 
into the proton component. We also assume that the gas is optically thin. 

The temperature of the gas determines the efficiency of the Coulomb energy transfer 
from the protons to the electrons and the rate of Bremsstrahlung cooling of the electrons. The 
balance between them defines whether the gas is in the two temperature regime (when the 
temperatures are high and the Coulomb collisions are very rare) or in the one-temperature 
regime (when temperatures are lower). 

In the two-temperature regime, we expect that the electrons are relativistic and the 
protons are non-relativistic: (?sel^ — kTe/inieC^) ^ 1, (?spl'^ = kTp/{mpC^) <^ 1. The 
viscous heating rate of the gas, the energy transfer rate from the protons to the electrons 
via Coulomb collisions, and the cooling rate of the electrons via Bremsstrahlung emission 
are given by 

gcoui = <5couiP^4^, Qcoui = 47rrglnA— ^, (8) 

Cse ^p 
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q~ = <5fr,R P^c^e, <5ff,R = 48a/rg— ^, (9) 

where af is the fine structure constant, Vf. is the classical electron radius, In A ^ 20 is the 
Coulomb logarithm, cj. ~ cj. , and we have neglected logarithmic corrections to the relativistic 
free-free emissivity. The subscript "R" in Qs^-r denotes relativistic Bremsstrahlung. 

In the one-temperature regime, both protons and electrons are cool and non-relativistic, 

and have nearly the same temperature, hence c^ ~ {me/mpjcl^ and c^ ^ 2c^ . The two 

energy equations (5) and (6) can be combined to yield a single energy equation for the 

accreting gas: 

pv del 2 dp + - r,n\ 

Cs^TK = (1-1 ' (10) 



-i-ldD ' dD 
where the free-free cooling takes the form 



1/2 2 

me c 



q = Qflf.NRP^C^, Qff.NR = StT ^^^afCTT ^ ^1^ , (11) 

m^ 

where gt is the Thompson cross-section, and the subscript NR stands for non-relativistic. 

As a result of high density of the gas in the BL, optically thin bremsstrahlung cooling 
dominates over self-absorbed synchrotron cooling; hence the latter may safely be neglected. 
We therefore neglect synchrotron emission in our analysis. In our model, we neglect the 
effects of radiation pressure compared to the gas pressure (which is of order M/MEdd *^ 1 
and, hence, negligible at low accretion rates) and Comptonization (which must be important 
in a high-temperature region, but is not important closer to the stellar surface, where the 
gas temperature is low, see more discussion below). For simplicity, we neglect also thermal 
conduction. Thus, our present model is similar to the models we used in our previous 
studies of hot accretion. Our simplified hydrodynamic model is, therefore, very instructive. 
It allows us to study the BL problem on the same grounds, on which the other hot fiow 
solutions have been treated, — as a viscous, radiative, purely hydrodynamic fiow. Once the 
basic hydrodynamic structure of the BL is understood, it will be worthwhile to put forward 
more sophisticated and detailed models, which should involve additional physics neglected 
in the present analysis. 

The set of equations (2)-(6) must satisfy certain boundary conditions at the neutron 
star. As the fiow approaches the surface of the star, the radial velocity must become very 
much smaller than the local free-fall velocity. For the fiow with f 7^ to match the radially 
non- moving stellar surface, the gas density should diverge, according to Eq. (2). Being 
proportional to the density squared, radiative cooling diverges as well, thus bringing the gas 
temperature to zero. In contrast, the angular velocity must approach the angular velocity 
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of the star Q^, hence remains nonzero. Naturally, the inner boundary conditions, as D = 
{R - /?,) -^ 0, are^ 

p-^oo, cl~^0, Vt^Vt^, V -^0. (12) 

The outer boundary conditions are set by the hot settling flow solution at i? ^ i?* (Medvedev 
& Narayan 2001). We treat the mass accretion rate M as a parameter. 

We would like to comment here that a self-similar solution is only an apptoximate 
solution far from the boundaries. In order to match boundary conditions, transition layers 
can form where variables relaxe from the values at the boundary to a self-similar behavior. 
These layers may be quite extended and be the site of interesting physical processes not 
accurately reproduced by the self-similar solution. As we shall see below, the self-similar 
solution for the hot boundary layer matches well the conditions at the inner boundary, 
whereas at the outer boundary the flow is required to settle onto a subsonic, hot inflow and 
a transition layer appears, see Figure 1 in the next Section. 



3. Self-similar solution for a boundary layer 

3.1. Two-temperature solution 

The gas in the two-temperature regime is governed by equations (2)-(6), which we now 
consider one by one and identify leading terms in them. 

Let us flrst consider equation (3). First, we note that the rotation of the gas is sub- 
Keplerian, fil <C ^"k^, so that we can neglect the flrst term in equation (4). Next, we cast 
the equation into the form: 

where fg^* = ^/2D.k^R^: is the free-fall velocity that near the stellar surface. We now make 
the following assumptions, which consistency with the obtained solution must be checked a 
posteriori: (i) the flow is always subsonic, f^ <^ c^, then the second term in the flrst brackets 
may be neglected and (ii) c^ grows with D slower than linearly (for c^ oc D^ we should have 



^These conditions correspond to an infinitely dense and cold star. To allow for a smooth match of a BL 
solution with a realistic, non-zero temperature star, an additional physics, such as thermal conduction, must 
be invoked. Such a consideration goes beyond the scope of the present study. 
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(3 < 1), then the second term in the second brackets is sub-dominant and may be neglected 
as well. With these assumptions, the equation simplifies to 

pel = p = const., (14) 

that is, the pressure is constant throughout the boundary layer. Note also that in the 
two-temperature regime, the gas pressure is dominated by the protons: 

2 I f^-Lp kle Tne\ kip 2 /t r-\ 

\ nip rrie rUpJ nip ^ 

because Tp^ T^, where Tp and Tg are the proton and electron temperatures. 

As we mentioned earlier, we are looking for the solution which is self-similar (i.e., power- 
law) in D. On the other hand, angular velocity Q approaches a constant at the star surface: 
r2 = f]^,. Thus, we readily conclude that fl oc D'^ (otherwise it is either zero or diverges at 
D = 0). Therefore, equation (4) ought to reduce to 

dn/dD = (16) 

with the right-hand-side of (4) being sub-dominant (i.e., it will introduce a small D-dependent 
correction to the zeroth order solution fi = f2*; we demonstrate this below). 

Since dVt/dD = 0, the heating rate in equation (5) vanishes. Together with the conti- 
nuity equation (2), the energy equations for the protons and the electrons read. 
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The system of equations (2), (14), (16)-(18) admits the following self-similar solution: 

p = Pod-'/\ c% = e%,d"\ cl = el,d'/\ v = vod^/\ n = nod', (19) 

where we used the dimensionless distance d = D/R^. The constant factors may be found as 
follows. Equation (17) determines the relation between c^q and po- The sum of equations 
(17) and (18), together with cKm^/nip) <^ c^p and Eqs. (19), yields the relation between 
c^o and po- On the other hand, at the outer edge, BL should smoothly match to the hot 
settling solution. This requires the pressure balance to hold at the interface, hence p = po = 
PqC^o = Pout- From this condition, all the constants are determined uniquely (the angular 



velocity is, of course, set by the stellar rotation): 



Po =P. 



'J:.{Q^,KQco.i/^y"'B'l\ (20) 



c 



= pI {Qs,KQconl/c') B-\ (21) 

do = PoQcoui5-\ (22) 

vo = Ml (47ri?2po) , (23) 

fio = fi*, (24) 



where we introduced 



M f 1 \ 

Finally, pout is estimated from the self-similar settling flow solution (Medvedev & Narayan 
2001) by setting the transition radius approximately equal to the radius of the star: 



16 \/Qs^YlQcou\Rg 

The transition from the boundary layer to the hot settling flow occurs at some distance 
dtri which can be estimated by matching the boundary layer density or sound speed (which 
is equivalent, because the pressure also matches) to that of the hot settling flow. Assuming 
that dir *^ 1, the hot settling flow has the proton sound speed squared is equal to (c^/6)r^^. 
Therefore, from c^g ~ (c^/6)r~^, we obtain: 

d,, ~ (6r,c^po/c^)-5/2 _ Q_Q44 ^27) 

for typical parameters, a = 0.1, rh = 0.01, s = 0.3, m = 1.4, r* = 3. 



3.2. One-temperature solution 

In the one-temperature regime, the temperatures of the protons and electrons are nearly 
equal, c^^ ~ {rrie/mpjcj.^, and the both species contribute equally to the pressure. The energy 
equations (17) and (18) reduce to the single energy equation for the accreting gas (10). 
Since dVt/dD = 0, the heating rate in equation (10) vanishes. Together with the continuity 
equation (2), the energy equation reads. 



M f I dci ci dp 



A-kRI \-i-l dD p dD 



Qs,nrP Cs {2t 
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The system of equations (2), (14), (16), (28), admits the following self-similar solution: 
p = po rf-2/5, d = clo d^/\ v = Vo rf2/5, Q = Qod°. (29) 

Here the constant factors are: 



. _ r.l/5n-2/5 r2/5 
yO — foutVff^NR-'-^ ' 

„2 _ „4^2 TD-2 



(30) 
(31) 



and Vq and i^o are given by equations (23), (24). 



3.3. Some comments 

Having derived the self-similar solutions, we now prove that the assumptions made in 
order to simplify equations (3) and (4) are consistent. 



Let us first consider the angular momentum equation (4). Because pc^ = Pout =const. in 
the boundary layer and angular momentum flux J =const. throughout the entire accretion 
flow, this equation has a simple solution: 



n 



n. 



fi* 



J 



MQM 



(1 - e-^/^) + e-^/^ 



J 



MQM 




n^{l + ed), 



where 



A = 4napout Rl/{nK,M) ~ 10-2a2g2^-i^i/2 ^ 
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(32) 



(33) 



and in the expansion we assumed D /A <^ 1. Note that for typical parameters of the hot 
settling flow: a ~ 0.1, s ~ 0.1, m < 0.01, r* ~ 3, we estimate A > lO'^^-R^,. To calculate 
the first-order correction to the solution, e, we use the expression for J from Eq. (15f) of 
Medvedev & Narayan (2001). Noting that J/{MVtMl) > 1, we obtain: 



(3/2)r;^ ~ (1/18) « 1; 



(34) 



note that e depends on r* only. Therefore the zeroth order solution describes the scaling of Vt 
accurately throughout the boundary layer (where D /A < 1), with the first order correction 
introducing just a minor effect. For progressively smaller NS angular velocities, f2*, the 
transition layer between the self-similar boundary layer solution and the hot settling flow, 
A < D < i?*, becomes more and more extended. Uniting the applicability of our self-similar 
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solution. Moreover, for small Q^, the expression for J used above becomes inaccurate, 
too. Hence we limit our investication to the case of a rapidly rotating neutron star. 

A similar analysis can be applied to equation (13). However, it is sufficient to prove 
that the omitted terms are sub-dominant. Then the ffist order corrections to the self-similar 
solution will also be small. Firstly, we readily see that the second term in the ffist brackets 
is small: v'^/c^ oc cP'^ ^ as rf ^ 0. Secondly, for the second term in the second brackets be 
sub-dominant, c^ should scale slower then linearly in d, i.e., if c^ oc d^ then (3 must be less 
then unity. From our solutions (19) and (29), we have (3 = 2/5 < 1. We should comment 
here that with dVt/dD ~ viscous heating is negligible in the boundary layer. From equation 
(28) it is clear that it is adiabatic compression of gas, dp/dD > that keeps the gas hot in 
the presence of Bremsstrahlung cooling. 

Finally, we conffim our theoretical results numerically. We use the numerical code 
used in our previous studies, which solves the system of the height-integrated viscous, hy- 
drodynamic equations using the relaxation method. The inhomogeneous grid was specially 
designed so that to resolve the thin boundary layer accurately. For more discussion, we refer 
the reader to our previous paper (Medvedev & Narayan 2001). The calculated structure of 
the boundary layer is presented in Figure 1. Note the remarkable agreement of this numeri- 
cal solution with the theoretical one: Q =const., p ?«const., and p, Tp, T^, v following the 
predicted scalings. 



4. Discussion 

It is instructive to summarize the scalings of various gas properties as functions of stellar 
and accretion parameters. 

p oc al/5s2/5^2/5^-l^-9/5^-2/5^ (35) 

Tp OC a4/5//5m-2/5mV;6/5rf2/5^ ^3g) 

V oc a-i/5s-2/5m3/5^0r;i/^rf2/5_ ^37) 

Here we remind that r* = R^/ Rg whereas d = D / R^ and we used that Tp oc c^^. 

The luminosity per logarithmic interval in D, i.e., dL ~ q~RldD ~ q~ RlD d{log D) , is 
then 

dL/d{\ogD) oc a^/^s^/^m^^/^mV^^/^rf^/s. (38) 

Hence, most of the energy is radiated far away from the NS surface. One may also note the 
strong dependence of dL/d(\ogD) on the NS spin and radius. The transition radius, Eq. 
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Fig. 1. — The structure of the boundary layer. Log of p (in 10 



18 pj^ 3^ 



T^^ Te (in Kelvin), 



r2/f2j^^,, f/c and p = pc^ (in dyne cm"^) as functions of Log{D / R^) are shown for rh = 
0.01, q; = 0.1, s = 0.3. Grayscale lines represent the self-similar solution. The temperatures 
and the velocity are sensitive to the value of the external pressure, which is estimated from the 
hot settling flow self-similar solution (Medvedev & Narayan 2001). The agreenemt becomes 
nearly perfect when pout is increased by a factor of 4, compared to Eq. (26). 
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(27), scales as 

d,, ~ (6r,c2po/c2)-5/2 oc a-^/^s-s/^m^/s^Ory^ (39) 

The total boundary layer luminosity is then estimated to be: 

Lbl oc a^'^s^'^rn^'^mh--^'^dli^ oc mmr::\ (40) 

which perfectly agrees with a simple energetic argument that Lbl ~ Mc^ /r ~ M-^^f^rh/r^, ~ 
mrhr"^ (Medvedev & Narayan 2001). 

Our self-similar solution is subsonic (otherwise a shock would form) with the Mach 
number of the infalling gas being 

M = v/cs oc rf^/^ (41) 

that is, if the gas is subsonic where the boundary layer matches the bulk accretion flow 
(which is always so, because the gas in the hot settling flow is highly subsonic), it will 
remain subsonic all the way down to the stellar surface. 

Next, we estimate the effect of Comptonization. The electron scattering optical depth is 
Tcs — pKcsD, where k^s = (Jt/it^p is the electron scattering opacity for ionized hydrogen. We 
can use our self-solution for p to calculate r^g. Alternatively, we may recall that the density 
matches that of the hot settling flow at the transition distance, and at in this region the 
electron scattering optical depth (Medvedev & Narayan 2001), Tcs — lO'^as^r"^ (assuming 
again that Dtr ^ R*), is < 1/3 for typical parameters a = 0.1, s = 0.1 and r* = 3. Using 
the scaling for the density, p oc (i~^/^, we estimate the optical depth in the boundary layer 
as 

res,BL ~ lO^'as^^'id/dtrf ~ (l/3)ao.i4i^;3 (^/^tr)'/' < 1, (42) 

where ao.i = a/O.l and similarly for other parameters. Hence, gas is optically thin to electron 
scattering for reasonable accretion parameters. The Compton ^/-parameter, y = 16(cse/c)^res, 
calculated for the boundary layer should match that of the hot settling flow at the transition 
distance. Using the result of Medvedev & Narayan (2001) and that c^g oc d^^^, we estimate 

1/ ~ 2 X 10^as\;\d/dtr) ~ 220ao.iSo.i^*":3(^Mr). (43) 

The effect of Comptonization is important in the high-temperature region of the boundary 
layer where y ^ 1, which occurs at distances larger or comparable to 

4 ~ 4.5 X 10--^«o:i%K3^tr. (44) 

Therefore, we conclude that Comptonization is not important deep inside the boundary layer, 
d < dc, and our one-temperature self-similar solution is accurate there. However, the electron 
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temperature profile in the two-temperature zone may be substantially modified by the inverse 
Compton scattering. Since, however, the electron-proton collisions are relatively rare (the 
plasma is two-temperature), other gas properties, such as the density, proton temperature, 
etc. should not be strongly affected. (Since the Coulomb collision rate increases with the 
decresing electron temperature, a more careful analysis is necessary. Such an analysis is 
beyond the scope this paper.) Note that in systems with slow rotation (s < 0.01) and 
with low-viscosity gas (a < 0.01), the structure of the boundary layer is not affected by 
Comptonization because dc > c?tr- 

Finally, we estimate the position of the photosphere, i.e., the distance at which the 
free-free optical depth becomes of order unity. The free-free opacity is approximately equal 
to ftff ~ 1.7 X 10^^^T~^/^(p/'mp)^ cm~^. Using our self-similar solution, we calculate the 
free-free optical depth: 

Tff = 5.0 X 10-25a-12/5s-24/5^11/5^-1^8/5^-6/5^ (45) 

The optical depth becomes greater than unity at distances below the photospheric radius: 

c^phot = 5.2 X 10-i«ao:?%t^aoi'^'^'<3' « 1- (46) 

Hence, the free-free emission in, essentially, the entire boundary layer is optically thin. 

5. Conclusions 

In this paper we presented the analytical self-similar solution describing the boundary 
layer, which forms in the vicinity of a spinning neutron star. Our solution is hot, highly 
subsonic and contains no shocks. The crucial difference of our work from others is that 
the bulk accretion flow is a quasi-spherical hot settling accretion flow, rather than a thin 
Shakura-Sunyaev disk. We believe, our solution may be relevant to a boundary layer of an 
advection-dominated accretion flow as well. We critically examined the limitations of our 
solution, especially the effect of Comptonization. We concluded that the one-temperature 
solution is accurate, whereas the two-temperature solution accurately represents the density, 
proton temperature, infall and angular velocity proflles, but the electron temperature proflle 
may be modifled by Comptonization. 

The author is grateful to Ramesh Narayan and an anonymous referee for insightful 
comments on the manuscript. This work has been supported by NASA and DoE. 
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